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FOREWORD 


It  is  of  continuing  interest  to  the  Navy  and  the  DoD  to  learn  something 
about  the  distribution  of  prime  numbers  and  also  to  explore  methods  of 
devising  prime-rich  equations.  This  report  is  a  continuation  of  work 
previously  done  with  a  view  to  contributing  to  the  use  of  prime  number  theory 
in  encription  and  other  fields. 

The  analysis  on  which  this  report  is  based  was  done  on  the  employee's  own 
time  and  it  is  being  published  by  the  Center  because  of  the  relevance  of  the 
subject  to  the  Navy. 
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INTRODUCTION 


1  7 

In  previous  reports1*4,  a  special  array  of  the  odd  positive 
integers  was -descri bed  for  which  solutions  of  the  Diophantine 
equation  I*x  -x+c  corresponded  to  the  columns  of  the  array  and 
solutions  of  I *xz+x-r  corresponded  to  its  rows;  c  and  r  are  each 

equal  to  2N-1  and  N*l,2,3 .  represents  both  the  column  and  row 

numbers  respectively.  An  infinite  set  of  “primitive  cell  arrays" 
( PC  A  *  s )  was  postulated  to  exist.  Each  such  array  is  derived  from 
the  special  array  and  is  limited  to  the  prime  values  of  c  (and  r) 
in  the  equations  above.  Table  1  shows  a  section  of  the  original 
array  with  prime  values  of  I  underlined  along  with  its  PCA  for 
c*r*3  shown  below  it.  It  was  proved1  that  the  PCA's  for  3  and  5 
are  infinite  in  extent  and  conjectured  that  they  exist  for  all 
prime  values  of  c  (and  r) .  Based  on  the  properties  of  these 
arrays  a  new  method  was  devised  for  finding  prime-rich  "column" 
equations  whose  densities,  Dc,  were  comparable  to  that  of  Euler's 
equation,  x  -x+41  for  which,  for  xaI  to  40,  all  values  of  I  are 
prime  i ,e  0^  =  100%.  The  method  of  analysis  involved  locating 
those  columns  for  which  no  value  of  I  was  =  0(mod  P),  P  is  prime, 
no  matter  how  far  the  column  extended  and  then  finding 
coincidences  of  columns  "empty"  of  3's  and  5's  i.e.  no  I's  =0(mod 
3,5)  and,  in  general,  columns  where  no  I's  are  congruent  to 

0(mod,3,5,7, . ,P).  Equations  and  columns  were  used 

interchangeably  as  are  rows  and  equations  in  this  paper.  Columns 
or  rows  are  called  prime-rich  if  the  density  of  primes  is  >  67.5% 
where  the  density  is  defined  as  Dc  *  100  x  (prime  values  of 
I)/(total  number  of  I's);  the  total  number  of  I’s  is  equivalent 
to  xN  and  will  be  considered  to  be  *  40  unless  otherwise 
stated.  The  columns  were  treated  at  some  length  elsewhere  ,  but 
the  rows  were  mentioned  only  briefly.  However  as  the  results  for 
the  rows  are  just  as  significant  and  interesting  as  are  those  for 
the  columns,  they  are  presented  here.  The  analysis  for  the  rows 
is  very  similar  to  that  developed  for  the  columns  and  will  not  be 
repeated  but  is  described  in  reference  2. 

ANALYSIS 

It  is  clear  that  for  the  columns  all  integers  (I's)  are 
positive  but  that  for  the  rows  some  values  of  the  I's  will  be 
negative.  This  means  that  for  some  prime-rich  equations  there 
may  be  duplications  of  primes;  for  example  for  I»xz+x-169  the 
negative  values  are  -167,-163,-157,-149,-139,-127,-113,-97,-79, 
-59,-37,-13  and  the  positive  values  are  13,41,71,103,137,173 
There  is  one  duplication  of  primes  i.e.  ±13.  The  density  0^9 
for  xn*40  is  90%  (36  out  of  40  values  are  prime)  If  the  negative 
value,  -13  is  excluded  is  87.5%.  However  as  x^  is  increased 

the  differences  in  the  two  percentages  diminishes  and,  for  x^ 
very  large,  becomes  negligible.  For  example  for  x^=2398  the  two 
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percentages  are:  40.8673...  if  the  duplicate  is  included  and 
40.8256...  if  it  is  excluded.  Another  instructive  example  is 
I*x^  +  x-7759  .  For  x*l  to  87  the  values  of  I  range  from  -7757  to 
-103;  for  x=88  to  124  the  corresponding  range  of  positive  values 
is  from  +73  to  +7741.  There  are  two  duplications  of  integers 
namely  +4987  and  +7003.  The  x  values  which  give  the  first  pair 
are  112  and  53  respectively  and  for  the  second  pair  121  and  27. 
However  4897  and  7003  are  composite  numbers!  This  means  that  all 
the  percentages  D7759  for  xN  =  40 , 100 ,200 , . . . , 10 ,000  and  beyond 
need  no  correction  with  respect  to  the  elimination  of  duplicate 
primes  no  matter  how  large  x»,  becomes.  The  correction  factor 
required,  if  any,  for  Table  4  will  be  discussed  below.  The 
computer  program  developed  for  determining  the  densities  of 
primes  is  the  same  as  that  described  previously  with  the 
exception  that  it  was  modified  to  record  the  absolute  values  of 
the  negative  I's. 

Table  2  shows  those  equations  which  describe  which  rows 
contain  no  integers  =0 ( mod  3);  which  have  none  =  Q(mod  5)  etc. 

For  example  rows  1,4,7,10, . 3N+1  contain  no  I's  =  0(mod  3). 

Rows  2,7,12, .  and  a  second  set  5,10,15,....  have  no 

I's  =0(mod  5).  By  solving  these  equations  e.g.  3Ni+l  and  5No+2 
and  also  3Ni+l  and  5^  +  5  those  rows  "empty"  of  I's  iOIrnod  3,5) 
can  be  found;  for  example  the  equations  (or  rows)  I=xz+x-(30N+13) 
and  I=x^+x- ( 30N+19 )  where  N=0 , 1 , 2 , 3 , . . .  are  for  for  rows  devoid 
of  I's  =  0(mod  3,5).  This  technique  can  be  extended  to  rows  for 
which  there  are  no  I's  s  0(mod  3 , 5 , 7 , 1 1 , . . . , P )  P  being  prime. 

Table  3  shows  the  result  of  carrying  out  this  procedure  to 
determine  all  rows  in  which  no  I's  are  =  0(mod  3,5,7,11).  All  of 
the  30  possible  combinations  are  shown  in  the  table  and  the 
results  of  the  solutions  for  each  combination  are  shown  in  the 
last  column.  For  a  comparison  with  a  similar  set  of  "column" 
equations  see  Table  1  of  reference  2. 

All  30  row  equations  were  evaluated  for  x  =  l  to  40  in  order 
to  compare  them  with  Euler's  equation  for  which  *100%.  Next 

N  was  varied  from  0  to  500  for  all-thirty.  The  program  which 
incorporated  a  "Fast  Prime  Search"3  program  as  a  subroutine  was 
designed  to  select  and  print  out  only  those  rows  out  of  the  15030 
examined  for  which  Dr  >  67.5%.  Eighty  five  equations  which  met 
this  criterion  were  found.  Only  16  of  these  are  listed  in  Table 
4.  These  are  the  ones  which  continued  to  show  high  Dr  values  as 
xN  was  increased  and  were  better  than  or  compared  favorably  with 
Euler's  equation  for  values  of  xN*2398  (where  ,°41  *50%  )  .  It  can 
be  seen  from  the  table  that  ten  of  these  equations  have  higher 
D  ' s  than  Euler's  for  x^>2398.  Figure  1  is  a  plot  for  a  few  of 
tnese  "curves"  of  Dr  (x  axis)  as  a  function  of  Xu.  Sixteen  of 
these  curves  were  averaged  and  the  results  are  also  plotted  in 
the  figure.  It  is  a  relatively  smoothly  varying  curve  as 
compared  to  the  individual  curves  shown. 

With  respect  to  the  fact  that  al 1  rows  but  the  first  contain 
some  negative  I's  it  turns  out  that  the  densities  Dp  for  all  but 
two  of  the  equations  of  Table  4  are  not  even  affected!  A  simple 
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two  of  the  equations  of  Table  4  are  not  even  affected!  A  simple 
program*  was  devised  to  compare  those  positive  and  negative 
values  of  I,  whose  absolute  values  were  identical  in  order  to 
determine  how  many  prime  pairs  existed  for  a  given  row.  From 
Table  4  it  can  be  seen  that  this  occurred  rarely.  The  last 
column  of  the  table  shows  how  many  duplicate  pairs  of  primes  did 
occur.  Of  the  sixteen  equations  only  two  were  found  to  have 
duplicate  primes  namely  r=6163  and  r=42,463.  The  former  had  one 
prime  pair  and  the  latter  had  one  prime  pair  and  one  composite 
pair  The  entries  Cl,  C2,  etc.  in  the  column  indicate  the  number 
of  composite  pairs  and  those  shown  as  PI  etc.  indicate  the  number 
of  prime  pairs.  This  means  that  no  corrections  for  duplications 
of  primes  need  be  made  for  any  but  the  two  equations  mentioned 
above  and  for  these  the  correction  is  negligible  for  large  xN  and 
is  easily  calculated.  In  column  x^=10,000  of  Table  4  the 
correction.if  prime  duplicates  are  excluded,  is  equal  to  -0.01%! 


CONCLUSIONS 

The  search  for  prime-rich  equations  of  the  type  I=x^-x+c  has 
been  extended  to  the  related  equations  I=x^+x-r  .  In  neither 
case  has  the  search  been  exhaustive  but  it  is  apparent  that  most 
such  equations,  up  to  c  and  r  values  of  "  1,500,000,  have  been 
found.  The  search  was  limited  to  those  equations  whose  densities 
were  >67.5.  Of  the  85  equations  which  met  this  criterion  most 
can  be  considered  prime-rich  but  fall  somewhat  short  of  the 
"standard"  chosen,  i.e.  Euler's  equation.  Nevertheless  they  do 
merit  future  study.  Similarly  there  are  quite  a  few  equations 
not  even  considered  which  have  densities  of  60  to  65%  which  may 
also  be  worth  examining.  One  example  is  I  =  x ^  +  x- 1  , 384 , 90 7  , 

0 r  =  65% .  Its  density,  for  x^  =  8000,  was  found  to  be  38.725%;  its 
projected  Df  for  xN  is  about  36%.  As  mentioned  in  a  previous 
report  anotner  facet  worth  examining  is  that  of  finding  equations 
which  are  very  prime-poor. 

An  additional  object  of  this  search  has  been  accomplished 
and  that  was  to  devise  a  method  of  determining  which  of  a  large 
number  of  equations  might  be  prime-rich.  It  was  done  by  finding, 
by  use  of  the  "Special"  array  and  its  derived  primitive  cell 
arrays,  those  columns  and  rows  ( cor respondi  ng  to  the  "c"  and  "r" 

equations)  which  contained  no  integers  =  Q(mod  3,5,7 . P) 

where  P  was  limited,  for  practical  reasons,  to  the  prime  values 
3,5,7,11. 

Finally  is  is  felt  that  the  work  described  above  has, 
perhaps,  contributed  a  bit  of  an  insight  into  the  somewhat 
organized  randomness  characteri Stic  of  the  distribution  of  prime 
numbers . 


★ 

see  Appendix  A 
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TABLE  1.  "SPECIAL"  ARRAY  Of  ODD  INTEGERS  WITH  PRIMES  UNDERLINED 
WITH  THE  DERIVED  PRIMITIVE  CELL  ARRAY  OF  3  SHOWN  BELOW  IT 


IC=X^-X+C 


C  =  1 

3 

5 

7 

9 

n 

13 

15 

17 

19 

xl  =  1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

r  - 

i 

1 

5 

11 

19 

29 

41 

55 

71 

89 

109 

3 

3 

9 

17 

27 

39 

53 

69 

87 

107 

129 

5 

T 

15 

TT 

37 

51 

TT 

85 

105 

ITT 

151 

7 

IT 

23 

35 

TT 

65 

TT 

103 

125 

TIT 

TTT 

9 

TT 

TT 

47 

63 

81 

ITT 

ITT 

147 

TTT 

201 

11 

31 

45 

TT 

79 

99 

TTT 

145 

171 

ITT 

229 

13 

TT 

59 

TT 

TT 

119 

143 

169 

197 

TTT 

TTT 

15 

57 

75 

95 

117 

141 

167 

195 

225 

257 

291 

17 

73 

93 

115 

139 

165 

ITT 

223 

255 

TTT 

325 

19 

TT 

113 

137 

16  3 

191 

T TT 

TTT 

287 

323 

361 

21 

111 

TTT 

it r 

ITT 

TTT 

251 

285 

321 

359 

399 

23 

133 

159 

187 

217 

249 

283 

319 

357 

397 

439 

25 

157 

185 

215  247  281  TIT 

Primitive  cell  array  for 

355 

c  = 

395 

p  =  3* 

TTT 

481 

_ 

_ 

3 

3 

- 

3 

3 

- 

3 

3 

- 

3 

- 

3 

- 

- 

3 

- 

- 

3 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

3 

3 

• 

3 

3 

- 

3 

3 

- 

3 

- 

3 

- 

- 

3 

- 

- 

3 

- 

- 

- 

- 

- 

- 

• 

- 

- 

- 

- 

_ 

3 

3 

- 

3 

3 

• 

3 

3 

• 

3 

- 

3 

- 

- 

3 

- 

- 

3 

- 

- 

- 

- 

- 

- 

• 

- 

- 

• 

• 

_ 

3 

3 

- 

3 

3 

- 

3 

3 

- 

3 

- 

3 

- 

- 

3 

- 

- 

3 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

★ 

I  n 

the  primi ti ve 

cell 

array 

(PCA) 

for 

every 

integer 

in  the 

"Special"  array  not  divisible  by  3  the  corresponding  position  in 
the  PCA  is  left  blank  whereas  each  integer  =  0(mod  3)  is  replaced 
by  the  numeral  3. 
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TABLE  2.  EQUATIONS  SHOWING  THE  NUMBER  OF  SETS  OF  EMPTY  ROWS  FOR 
EACH  PRIME  VALUE  OF  P  FOR  THE  PRIMITIVE  CELL  ARRAYS 


p 

No.  Of 

Sets  of  "empty"  rows 

sets 

per  array  (N=0,l,2,..) 

3 

1 

3N+1  * 

5 

2 

5N+2  ,•  5N+5 

7 

3 

7N+1,  7N+2,  7N+6 

11 

5 

11N+2,  11N+3,  11N+4 ,  11N+8,  11N+11 

13 

6 

13N+2,  13N+3,  11N+5,  13N+6,  13N+11,  13N+12 

17 

8 

17N+1 ,  17N+4,  17N+5,  17N+6,  17N+8.17N+14.17N+16  17N+17 

19 

9 

19N+2 ,3,4,5 ,7 ,9 ,14 ,15 ,18 

23 

11 

23N+1, 3, 5, 6, 7, 8, 14, 16, 19, 20, 23 

29 

14 

29N+3, 5, 8, 9, 10, 11, 13, 17, 20, 24, 26, 27, 28, 29 

31 

15 

31N+2, 3, 4, 5, 7, 8, 10, 11, 15, 18, 23, 24, 25, 27, 29 

37 

18 

37N+1, 2, 4, 6, 8, 9, 12, 14, 15, 16, 17, 21, 26, 30, 31, 32, 33, 35 

p 

(P-D/2 

PN+M1,M2,M3 . ’M(P-l)/2 

*Subscripts  such  as  1,2,3,  etc.  in  3N3+1,  5N2+2, .  as  shown  in  Table  3  and 

elswhere  were  omitted  for  the  sake  of  clarity 


!j 
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TABLE 

3. 

COINCIDENCES 

OF  EMPTY 

ROWS  FOR  DIVISORS  3  AND 

5;  AND  3 

,5,7;  AND 

3,5,7,11  FOR  DERIVED 

PRIMITIVE 

CELL  ARRAYS 

No. 

of 

Equations  for 

Fami  I ies 

of  equations 

K  family 

"empty" 

rows  for: 

of  form 

x^+x-r:  r= 

1 

3!)i+1 

5N2+2 

7N->+l 

ll  ^ 

llN4+2 
"  3 

2310N+  883 

2 

ii 

11  1303 

3 

II 

ll 

ll 

"  4 

ll 

"  1723 

4 

II 

II 

II 

"  8 

ll 

"  1093 

5 

II 

ll 

II 

"  11 

ll 

"  43 

6 

II 

ll 

7N^+2 

11N.+2 

II 

"  1543 

7 

II 

ll 

1,4  3 

II 

"  1963 

8 

II 

ll 

it 

"  4 

ll 

"  73 

9 

II 

It 

a 

"  8 

ll 

"  1753 

10 

II 

ll 

ii 

"  11 

ll 

"  703 

11 

II 

ll 

7No+6 

llN4+2 

ll 

"  1873 

12 

II 

ll 

n J 

"  3 

M 

"  2293 

13 

II 

ll 

ii 

"  4 

H 

"  403 

14 

ll 

H 

M 

"  8 

ll 

"  2083 

15 

II 

II 

n 

"  11 

ll 

"  1033 

16 

II 

5N2+5 

7N^+1 

11N.+2 

II 

"  2269 

17 

11 

n° 

"  3 

ll 

"  379 

18 

II 

ll 

ii 

"  4 

ll 

"  799 

19 

II 

ll 

ii 

"  8 

ll 

"  169 

20 

II 

II 

ii 

"  11 

II 

"  !  42  9 

21 

II 

ll 

7No+2 

ii J 

llN4+2 

it 

"  619 

22 

II 

II 

"  3 

>1 

"  1039 

23 

II 

II 

ii 

"  4 

ll 

"  1459 

24 

II 

II 

ii 

"  8 

ll 

"  829 

25 

II 

II 

ii 

"  11 

ll 

"  2089 

26 

II 

II 

7N-J+6 

1  IN/, +2 

ll 

"  949 

TABLE  4.  EQUATIONS  x2+x-r  FOR  r=501229,  etc.  WHERE  THE  DENSITY,  D  IN 
PERCENT  (100xPRIMES/40)  FOR  x=l  TO  40  IS  >  67.5%  AND  WHERE  D_ 
IS  ALSO  GIVEN  FOR  xN  =100,  200,  ETC. 


Equations  having  higher  Dp  values  then  c=41*  for  xN=  10,000 


xni40 

100 

200 

400 

1000 

2398  5000 

8000 

10000 

±Pai rs 

501229 

67.5 

62 

63.5 

62.5 

67.5 

60.34..  54.3 

50.625 

49.250 

Cl 

349513 

67.5 

68 

65 

65.75 

68 

59.96..  53.62 

50.275 

48.85 

0 

445473 

67.5 

63 

59.5 

62.5 

65.9 

58.84..  52.68 

49.65 

47.95 

Cl 

800629 

67.5 

56 

59.5 

58.25 

61.1 

56.96..  51.16 

47.8375 

46.87 

0 

249439 

72.5 

67 

64 

64.75 

65.5 

56.79..  51.06 

48.1 

46.65 

0 

82009 

72.5 

65 

64.5 

68.75 

60 

51.41..  46.8 

44.475 

43.32 

Cl 

90073 

70 

67 

67 

71 

61.5 

53.79..  48.82 

44.8875 

44.7 

Cl 

152839 

75 

59 

61.5 

65.5 

60.8 

52.87..  47.74. 

45.05 

43.71 

0 

366349 

67.5 

56 

52.5 

54.25 

57.3 

51.79..  47.0 

43.625 

42.19 

0 

1063903 

67.5 

50 

46 

49 

50.7 

50.25..  45.84 

43.325 

42.1 

C2 

41* 

100 

86 

78 

67.5 

58.1 

50  45.22 

45.575 

41.49 

- 

Equations 

with  Dr 

values 

close  to  those 

for  r=41 

53509 

75 

63 

66.5 

68.5 

56.3 

49.58..  43.68. 

41.325 

40.11 

C2 

7759 

67.5 

77 

70.5 

62.5 

54.5 

47.83..  43.18. 

40.525 

39.25 

C2 

2293 

77.5 

77 

68.5 

63 

54.5 

47.28..  42.94. 

40.3625 

39.25 

PI 

6163 

70 

75 

67.5 

61 

53.8 

48.12..  43.1 

40.375 

39.16 

Cl, PI 

23083 

67.5 

67 

71.5 

63.75 

54.3 

47.70..  42.8 

40.375 

39.15 

PI 

42463 

70 

60 

63 

61.5 

54 

47.83..  42.68. 

39.95 

38.94 

C2 

*Euler's 

equation, 

x2-x+41,  a  "column" 

equation,  is  included  here  for 

the  sake 

of  comparison  with  the  row  equations. 
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Appendix  A 

If  only  odd  integer  solutions  of  the  equation  I=x^+x-r  are  considered 
then  it  is  obvious  that  x  and  r  must  be  integers.  If  there  exists  a  pair  1^, 
U  such  that  Ij»-I2  then  2  equations  exist  i.e.  Ij^x^+xj-r  and  I2=x22+x-r.i 
Adding  the  two  together  yields: 

1.  x^+x^+(x2^+x2~2r)=0 

and  solving  for  x,, 

2.  x1»l/2(-l+  /{8r+l-4x2[x2+l]  }) 

where  the  minus  value  of  the  square  root  sign  is  ignored  as  all  x^,  x2  values 
are  positive.  1 

For  a  pair  of  integers  x^,  x2  to  exist  such  that  I^=- Ig  the  right  hand 
side  of  equation  2.  must  be  an  integer  i.e.  the  expression  in  parentheses  is 
even  and  that  within  the  square  root  sign  is  an  odd  square.  If  x^  is  not  an 
integer  for  any  value  between  1  and  /r  then  there  are  no  duplicate  I's  whether 
prime  or  composite. 

Another  program  beside  the  one  shown  below  was  used  to  search  for 
duplicate  I's.  It  involved  examining  only  those  values  of  x ^  between  /r 
and  /(2r)  as  there  are  fewer  positive  values  of  I  in  this  range  than  in  the 
corresponding  range  of  negative  I's  in  the  range  x=l  to  /r.  Both  programs 
were  designed  to  indicate  and  print  out  only  those  matching  pairs  of  I's  (and 
pertinent  x  values)  which  were  found.  The  somewhat  simpler  but  slighly  longer 
program  shown  below  was  used  to  confirm  the  results  of  the  one  mentioned  above 

The  Basic  program  run  on  tne  HP  983 6C 


10 

!  This  is  primedupl  1,  of 

3/9/87 

20 

OPTION  BASE  1 

30 

DIM  Kr( 16) 

40 

DATA  2293,6163, . 

,10639 

50 

READ  Kr(*) 

60 

FOR  Rr=l  TO  16  STEP  1 

70 

Rp=Kr(Rr) 

80 

Sr=SQR( Rp) 

90 

Si*INT(Sr) 

100 

IF  (Sr-Si )>.5  THEN 

110 

Sr=Sr+l 

120 

END  IF 

130 

FOR  X*1  TO  Sr- 1  STEP  1 

140 

P=SQR(  8*Rp+l-4*X*( X+l) ) 

150 

Pi  =  I  NT(  P) 

160 

Ps*P-Pi 

170 

IF  Ps>0  THEN  GOTO  200 

180 

IF  Ps*0  then 

190 

PRINT  "P="  ;P;"X="  ;X;"THEKE 

IS  A 

191 

ENO  IF 

200 

NEXT  X 

201 

NEXT  Rr 

210 

END 

■RP=" ;RP 


*  See  the  first  column  of  Table  3  (value  41  not  included  in  40  above) 
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